VARIABLE STARS IN GALACTIC GLOBULAR CLUSTERS 
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ABSTRACT 

Based on a search of the hterature up to May 2001, the number of known 
variable stars in Galactic globular clusters is approximately 3000. Of these, 
more than 2200 have known periods and the majority (approximately 1800) 
are of the RR Lyrae type. In addition to the RR Lyrac population, there 
are approximately 100 eclipsing binaries, 120 SX Phe variables, 60 Cepheids 
(including population II Cepheids, anomalous Cepheids and RV Tauri) and 120 
SR/red variables. The mean period of the fundamental mode RR Lyrae variables 
is 04585, for the overtone variables it is 04342 (04349 for the first-overtone 
pulsators and 0^296 for the second-overtone pulsators) and approximately 30% 
are overtone pulsators. These numbers indicate that about 65% of RR Lyrae 
variables in Galactic globular clusters belong to Oosterhoff type I systems. 
The mean period of the RR Lyrae variables in the Oosterhoff type I clusters 
seems to be correlated with metal abundance in the sense that the periods are 
longer in the more metal poor clusters. Such a correlation does not exist for 
the Oosterhoff type II clusters. Most of the Cepheids are in clusters with blue 
horizontal branches. 

Subject headings: stars: horizontal-branch — stars: oscillations — stars: 
variables: RR Lyrae, Cepheids, SX Phe — binaries: eclipsing — globular 
clusters: general 
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1. INTRODUCTION 

Variable stars are useful standard candles for determining the distances to the Galactic 
globular clusters (GGCs) and it is important to know these distances for an understanding 
of the age, structure and formation of the Galaxy. As a result, there are many papers on 
the subject of variable stars in globular clusters. More than sixty years ago. Sawyer (1939) 
pubhshed a catalog of variable stars in GGCs to enable investigators interested in the topic 
to get a clear picture of exactly what work had been done in the area. The catalog listed 
1116 variables in 60 clusters. Second and third editions of the catalog, with 1421 and 2119 
entries respectively, were subsequently produced (Sawyer 1955, Sawyer Hogg 1973). Sawyer 
Hogg intended to publish a fourth edition, and in preparation for this, she recorded material 
from the relevant papers published between 1973 and 1988 on reference cards. After her 
death in 1993, the cards were placed with the University of Toronto Archives and Records 
Management Services (UTARMS) where they now reside. In 1997, an electronic version of 
the 1973 catalog was produced and the material from the 1973-1988 reference cards was 
included (Clement 1997). 

Since 1988, there have been many papers on the subject of globular clusters and 
their variables. Furthermore, during this period CCD detectors have been widely used. 
Consequently, many new variables have been discovered, particularly in crowded central 
regions of clusters. Thus it is an appropriate time to update the catalog and prepare a 
summary of the material. To search the literature for papers published after 1988, we 
consulted volumes 49-58 and 61-68 of the Astronomy and Astrophysics Abstracts for 
the years 1989-1993 and 1994-1997 respectively, and for papers published after 1997, we 
consulted the NASA Astrophysics Data System. In this paper, we summarize the numbers 
and types of variable stars in GGCs. In addition, for the RR Lyrae variables and Cepheids, 
we illustrate how the periods relate to the clusters' metal abundance and horizontal branch 
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morphology. The complete updated catalog can be obtained at the following website: 
http:/ /www. astro. utoronto.ca/people.htmi 

2. SUMMARY OF THE DATA 

In Table 1, we list the Galactic globular clusters which are known to contain 
variables. For each cluster, we list the horizontal branch ratio (HBR) and the metal 
abundance, [Fe/H]. The HB ratio, {B — R)/{B + V + R), is a parameter devised by 
Lee (1990) to describe horizontal branch morphology. Both of these quantities were 
taken from the 1999 update to Harris's (1996) catalog of globular cluster parameters 
(http : //www.physun.physics.mcmaster.ca/Globular.html). The total number of variables 
in each cluster and the number for which periods have been determined is listed in column 
(3) and the main types of variables are indicated in columns (4) to (9). These totals do 
not include any stars that are suspected or confirmed to be field stars. 'SX' (column 4) 
refers to the SX Phe type variables, stars in the region just above the ZAMS, that pulsate 
with periods less than 0^2. 'RR' (column 5) refers to RR Lyrae variables. Although many 
variables have been tentatively classified as RR Lyrae, we include here only the ones for 
which periods have been determined. 'Cep' and 'RV Tau' (column 6) refer to Cepheids 
(anomalous or population II) and RV Tauri variables. In general, we have classified stars 
with periods less than 1 day as RR Lyrae variables and stars with periods greater than 1 day 
as Cepheids. However, the actual threshold between these two groups of variables may be 
P=0'!75 or O'^.SO (Wallerstein & Cox 1984, hereafter WC84; Gautschy & Saio 1996, hereafter 
GS96). The 'SR/Red' variables (column 7) are pulsating variables that have periods greater 
than 35 days and/or stars that have been classified as semi-regular or irregular. Column 
(8) lists the number of eclipsing binaries and in column (9), we cite the papers in which 
new information has been presented concerning the elements of the variables since Sawyer 
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Hogg's (1973) catalog. For most of the clusters, the totals of columns 4 to 8 are equal to or 
greater than the number of variables for which periods have been determined because many 
irregular variables do not have published periods. However, in some cases (e. g. M15), there 
are 'unclassified' variables with published periods. 

There are 147 clusters in the Harris catalog, but only 102 are known to have variable 
stars. This does not mean that there are none in the other 45 because not all clusters 
have been searched. Also, with the use of CCD detectors and new search techniques, for 
example, the image subtraction method introduced by Alard & Lupton (1998), many new 
variables will be discovered, even in the clusters that have already been studied. 

In Figure 1, wc plot the period-frequency distributions for pulsating variables in 
globular clusters. In the top panel, stars with periods less than 1 day are plotted. The 
stars with P<0l2 are considered to be SX Phe variables. Most of the others are RR Lyrae 
variables, but a few of the longer period ones may be anomalous Cepheids or population II 
Cepheids. The middle panel includes the Cepheid and RV Tauri variables with periods in 
the range 1 to 50 days. The RV Tauri stars are plotted with their pulsation periods. The 
plot also includes two stars (one in cu Centauri and one in NGC 5466) with periods less than 
1 day and classified as anomalous Cepheids. In the lower panel, variables with pulsation 
periods between 35 and 440 days are shown. 

Figure 2 shows the period-frequency distribution for the eclipsing binaries. The 
majority of these stars have periods less than 1 day and about half are of the W Ursa 
Majoris type. Rucinski (2000) recently reviewed the hterature on W UMa-type binaries 
and found that at least one-third of his sample were foreground stars. We have not included 
any of these foreground stars in Figure 2 or in Table 1. 
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3. THE SX PHE VARIABLES 

Rodriguez & Lopez- Gonzalez (2000) published a catalog of SX Phe stars in GGCs and 
in nearby galaxies. They showed that there seems to be a correlation between mean periods 
and metallicity in the sense that the periods of the SX Phe variables are longer in stellar 
systems with higher metaUicity. They also pointed out that all of the SX Phe variables 
have been discovered in the last 20 years and that there has been a great increase in their 
numbers during the last few years. As more clusters arc studied, the number of SX Phe 
stars will undoubtedly increase further and we will be able to gain a better understanding 
of how their properties relate to the properties of the clusters to which they belong. An 
interesting point to note from the data of Table 1 is that, in four of the 21 clusters known 
to have SX Phe stars, no RR Lyrae stars have been detected. Three of these, NGC 4372, 
6397 and 6752, have extremely blue horizontal branches and the other one, M71, has an 
extremely red HB. Since SX Phe stars are near the main sequence, it seems that they are 
not affected by horizontal branch morphology. 



4. THE RR LYRAE VARIABLES 

For the RR Lyrae variables, we have adopted a system of notation recently introduced 
by the MACHO consortium (Alcock et al. 2000, hereafter AOO). 'RRO' instead of 'RRab' 
is used to designate fundamental mode pulsation, 'RRl' instead of 'RRc' for first-overtone, 
'RROl' instead of 'RRd' for double-mode (fundamental and first-overtone), 'RR2' instead of 
'RRe' for second-overtone pulsation and 'RR12' for double- mode (first and second-overtone). 
No RR12 variables have yet been identified in globular clusters, but AOO found three in 
the LMC. In Table 2, we list information concerning the periods of RR Lyrae variables in 
the GGCs: in columns (2) and (3), the number of fundamental mode pulsators and their 
mean periods, in columns (4) and (5), the number of first overtone pulsators and their 
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mean periods, in column (6) and (7), the number of second overtone pulsators and their 
mean periods, in columns (8) and (9), the number of double-mode pulsators and their 
mean periods and in column (10), the mean fundament alized period. Since the dominant 
mode of pulsation for most of the double-mode variables is the first-overtone, we list the 
mean overtone periods in column (9). The fundament alized period was computed assuming 
period ratios: Pi/ Pq = 0.745 and P2/P1 = 0.804, typical period ratios derived by AOO for 
the RROl and RR12 variables in the LMC The possibilty that RR2 variables really exist 
has been the subject of some discussion in the literature. In recent studies of u Cen and 
NGC 5897, Clement & Rowe (2000, 2001) showed that the RRl and RR2 variables seem to 
separate into two sequences in the period-luminosity and period-amplitude plots. Therefore, 
to decide whether an overtone pulsator should be classified as RRl or RR2, we consulted 
the period-amplitude relation of Clement & Rowe (2000). However, it is not always easy 
to discriminate between RRl and RR2 variables in the period-amplitude plot because the 
models of Bono et al. (1997) predict that the period-amplitude relation for RRl stars shows 
a characteristic 'bell' shape with amplitudes decreasing at shorter periods. Thus in some 
cases, our classifications are uncertain. 

According to Table 2, the mean period for all of the RRO variables in GGCs is 04585 
and for the overtone variables, it is 04342 (04349 for the first-overtone and 04296 for the 
second-overtone). Approximately 30% of the RR Lyrae variables are overtone pulsators. In 
calculating the mean overtone periods, we have included the RROl variables as first-overtone 
pulsators because the first-overtone is generally their dominant mode of pulsation. These 
figures indicate that most of the GGC RR Lyrae variables are in Oosterhoff type I systems.0 

^ Oosterhoff (1939, 1944) recognized that globular clusters could be classified into two 
groups according to the period-frequency distribution of their RR Lyrae variables. He showed 
that the mean periods for the RRO variables were ~0455 in the type I (Ool) clusters and 
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In their investigation of the MACHO data for LMC RR Lyrae variables, Alcock et ah 
(1996, hereafter A96) reached a similar conclusion. 

In Figure 3, we show the period- frequency distributions for the RR Lyrae variables 
separated according to Oosterhoff type. In assigning the Oosterhoff type, we assume that 
all clusters for which the mean period of the RRO stars is <0460 belong to group I and 
that clusters for which the mean RRO period is >0460, with the exception of Ruprecht 
106, belong to group II. The study of the RR Lyrae variables in Rup 106 by Kaluzny et 
al. (1995) shows that its period-amplitude relation is similar to that of M3 and so we have 
classified it as Ool. Pal 5, which has no RRO stars, is considered to be Ool on the basis 
of the mean period of its RRl variables. When the period-frequency data are separated 
by Oosterhoff type, it is possible to see that there are two peaks in the distribution for 
the overtone pulsators, particularly in the OoII clusters. This is an effect that A96 noted 
in the MACHO LMC data for RR Lyrae variables and they attributed it to pulsation in 
the second-overtone mode. The mean period for the 936 RRO stars in the Ool clusters is 
0'!559, for the 232 RRl and RROl stars, it is 0'!326, for the 29 RR2 stars, it is 0'!281, the 
mean fundament alized period is 04533, and 22% of the variables are overtone pulsators. 
The equivalent figures for the OoII clusters are 0^659 for 333 RRO, 0^368 for 268 RRl and 
RROl, 0^306 for 44 RR2, 0^580 for < Pf > and 48% are overtone pulsators. These numbers 
are comparable to what Oosterhoff (1939) found for the two groups, but the percentage of 
the (lower amplitude) overtone pulsatorsQ has increased. 

~0'!65 in the OoII clusters. He also found that, in the OoII clusters, the percentage of 
overtone pulsators was higher than in the Ool clusters. 

^When Oosterhoff discovered the two cluster groups in 1939, he used the terms 'c-type' 
and 'a- and b-type' for first-overtone and fundamental-mode pulsators. This is because it 
was not recognized until the following year that the c-type variables were pulsating in the 
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In Figures 4 to 7, we plot [Fe/H] and HBR versus the mean periods of the RRO and 
RR10 variables. Since any correlations are more significant if the sample of RR Lyrae 
variables is greater, we plot in the upper panels only the systems with at least 15 RRO stars 
(Figures 4 and 5) or at least 15 RRl stars (Figures 6 and 7) and then in the lower panels we 
include systems with fewer variables. Similar plots for the mean 'fundamentalized' periods 
are shown in Figures 8 and 9. In Figure 4, there seems to be a correlation between [Fe/H] 
and mean period for the Ool clusters (i.e. the ones with < PO > less than 0460) in the 
sense that the periods are longer in the more metal poor clusters. However, the correlation 
breaks down when the mean period is greater than 0460, i.e. in the OoII clusters. A similar 
trend can be seen for the mean fundamentalized periods in Figure 8. The mean period is 
correlated with [Fe/H] for clusters with < Pf > less than about 0456 (Ool clusters), but not 
for clusters with longer periods. The period-[Fe/H] correlation can also be seen among the 
RRl mean periods less than 0434 in Figure 6, but it is not as marked as in Figures 4 and 8 
because Ool clusters do not have so many RRl variables. 

Until recently, it has been assumed that metal rich clusters have RR Lyrae variables 
with short periods and therefore belong to the Ool class, but investigations of NGC 6388 
and NGC 6441 by Pritzl et al. (2000) and Layden et al. (1999) have shown that these two 
clusters are exceptions to the rule. Both have [Fe/H]~ —0.6 which makes them the most 
metal rich clusters with RR Lyrae variables, but their < PO > values are greater than 0470. 

first-overtone mode (Schwarzschild 1940). 

^No plots were made for the RR2 stars because there are not enough of them, u Cen is 
the only cluster in which we have classified more than 5 stars as RR2. We considered the 
possibility that our assumed mean periods for the RRl stars would be in error if the stars 
we considered to be RR2 variables were in fact RRl variables. However, it turned out that 
such errors would not be large enough to affect the appearance of Figures 6 and 7. 
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Because of these long periods, we classify these clusters as OoII. The discrepant 'crosses' in 
the middle and lower panels of Figure 4, the lower panel of Figure 6 and the middle and 
lower panels of Figure 8 are the points for these clusters. In fact, NGC 6441 does not even 
appear in Figure 8. Since it has no known RRl stars, its mean fundamentalized period is 
so long (04768) that it is off the scale of the diagram. The metal rich cluster 47 Tuc (NGC 
104) also has an RRO variable with a period longer than 0^70 (Carney et al. 1993), but it 
does not appear in Figures 4 or 8 because it is the only RR Lyrae variable in the cluster. 

The importance of [Fe/H] for determining properties of globular cluster RR Lyrae 
variables has already been questioned by Clement & Shelton (1999, hereafter CS99) and 
by Lee & Carney (1999b). Lee & Carney compared M2 and M3, two clusters with similar 
metal abundance and found substantial differences in their RR Lyrae populations. CS99 
studied the period ^-amplitude relation for RRO stars with 'normal' light curve^ in several 
globular clusters and found that the V amplitude for a given period was a function of 
Oosterhoff type, but not a function of [Fe/H]. The OoII clusters in their sample were M9 
and M68, but they also included two M92 RRO stars for which published CCD photometry 
was available. Subsequently, Clement (2000a) showed that the RRO variables in in the 
OoII cluster M55 and in the metal rich clusters 47 Tuc and NGC 6441, as well as the 
RRO stars brighter than V = 14.65 in uj Cen (the OoII variables) all fit the OoII relation 
of CS99. In the meantime, Papadakis et al. (2000) confirmed that the RRO variables in 
the OoII cluster NGC 6426 fit the OoII period-amplitude relation derived by CS99. More 
recently, Kopacki (2001) compared the period-amplitude relations for the RRO variables in 
three OoII clusters with different [Fe/H], M2, M53 and M92, and verified that there was no 
significant period-shift. Thus period-shift for the OoII clusters is not a function of metal 

^The compatibility test of Jurcsik & Kovacs (1996) was used to assess whether or not the 
light curve of a star was normal. 
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abundance. The P-A relation plotted by Alves et al. (2001) for NGC 5986 also seems to 
coincide with the CS99 OoII relation. An interesting feature of their diagram is that the 
P-A relation for NGC 5986 is shifted to longer periods than the fiducual ridge line they plot 
for M15. Since M15 is more metal poor than NGC 5986, they were expecting the M15 RRO 
stars to have longer periods. M15 which is often considered the prototype for OoII clusters 
seems to have a P-A relation that is different from the others! As noted above, Figures 4, 6 
and 8 also illustrate that the periods for OoII clusters do not depend on [Fe/H]. It therefore 
seems well estabhshed that, for the OoII clusters, mean periods and period shifts are not 
correlated with [Fe/H]. 

However, the situation is different in the Ool clusters. We have seen from Figures 4 
and 8 that there is a correlation between mean period and [Fe/H] for the Ool clusters. The 
Ool clusters that CS99 used for their study were M3 amd M107. They acknowledged that 
the transition between fundamental and first-overtone mode pulsation occurs at a shorter 
period in M107 compared to M3. However, the short period fundamental mode pulsators in 
M107 did not have 'normal' light curves. As a result, there was no period shift between the 
period-amplitude relations for 'normal' RRO stars in the two clusters. Later, Kaluzny et al. 
(2000) performed a similar analysis for the RRO variables in M5 and found that its P-A 
relation was shifted to shorter periods compared with M3. Borissova et al. (2001) reached 
the same conclusion in a recent study of NGC 6229. Thus in the Ool clusters, mean periods 
and period shifts seem to be correlated with [Fe/H], but in the OoII clusters they are not. 
This is another important difference between clusters of the two Oosterhoff groups and may 
indicate something about the evolutionary status of their RR Lyrae variables. 

In Figures 5, 7 and 9, HBR is plotted against mean period. In Figure 5, we see 
that clusters with [Fe/H]< —1.6 and blue horizontal branches (HBR> 0) belong to the 
Oosterhoff type II class. However, having a blue HB does not necessarily guarantee an 
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OoII classification because there are Ool clusters with HBR> 0.5. On the other hand, no 
cluster with HBR< belongs to the OoII class.0 According to Lee et al. (1990), the RR 
Lyrae variables in OoII clusters have evolved away from the ZAHB, and as a result, have 
longer periods and higher luminosities than the ones in Ool clusters, which are ZAHB stars. 
Furthermore, Lee (1990) calculated models illustrating that for metal poor clusters with 
HBR> 0.65, the bluer the HB morphology of the cluster, the longer the periods should be. 
To test this hypothesis, Catelan (1994) made plots of mean periods and HBR using the 
cluster data available at the time and found no correlation. It is interesting to note that 
with the new data plotted in Figure 5, the cluster with the longest mean RRO period (M2) 
is the one with the bluest horizontal branch. Another feature to note in these diagrams, 
particularly Figure 9, is that when all the clusters are considered, there is a weak correlation 
between mean period and HBR. 

5. THE POPULATION II CEPHEIDS 

The period-frequency distrubution for the Cepheids and RV Tauri variables is shown 
in the middle panel of Figure 1. The four RV Tauri variables, plotted with their pulsation 
periods which range from about 29 days for VI in uj Cen to 46 days for VI 7 in M28 seem to 
form a long period extension to the Cepheids. The gap in the period at around 10 days has 
an evolutionary origin. The longer period variables are either on a blueward excursion from 
the AGB during He-shell flashes or they are on the way towards the white dwarf cooling 
region (GS96, Gingold 1976). Their periods range from about 12 days (V17 in M14) to 29 

^The cluster with HBR=-0.8 and < P{RRO) >, < Pf >= 0.617 is Rup 106 which we 
have classified as Ool. This anomalous cluster was shown to be extremely young by Chaboyer 
et al. (1992). 
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days (V21 in M28). 

The stars with periods less than 10 days have often been referred to as BL Her 
stars, but GS96 pointed out that this is not an appropriate name because BL Her itself 
is not metal deficient. They prefer to call them AHBl stars because they have evolved 
from the blue HB and are on the way to the AGB. They cross the instabilty strip at 
higher luminosities than the RR Lyrae variables and therefore have longer periods. It is 
expected that their periods should increase with time and Wehlau & Bohlender (1982) have 
confirmed this to be the case. Four of the stars in the short period group are considered to 
be anomalous Cepheids. These are V19 in NGC 5466 (Zinn & Dahn 1976) and three uj Cen 
stars classified by Kaluzny et al. (1997c). According to WC84, anomalous Cepheids (AGs) 
have periods between about O'ISO and 3 days, but are too luminous for their periods. The 
high luminosities can be accounted for if they have larger masses than other population 
II Gepheids and these large masses are thought to occur because they are formed by the 
coalescence of one or more stars. Two of the globular cluster AGs in our sample (Ogle # 
161 in u) Cen and V19 in NGC 5466) have periods less than a day. 

In Figure 10, we plot [Fe/H] and HBR versus period for the individual Gepheids and RV 
Tauri variables. There is no correlation between period and [Fe/H], but most have [Fe/H] 
between —1.25 and —1.75. In fact, all but one of the Gepheids are found in clusters that 
are more metal poor than —1.25. Horizontal branch morphology is an important parameter 
for determining whether or not a cluster will have Gepheids. Generally they occur only in 
clusters with blue horizontal branches and this can be seen in Figure 10. However, we also 
see in Figure 10 that two clusters with HBR< have Gepheids. These are NGC 2808 and 
Pal 3, and in both cases the Gepheids have periods less than 4 days. Another interesting 
feature of the diagram is that, among the Cepheids in clusters with HBR less than 0.5, the 
periods tend to be longer in clusters with bluer horizontal branches. 
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It is often assumed that the boundary between RR Lyrae variables and Cepheids occurs 
at a period of one day. However as we noted in section 2, it is probably shorter. WC84 
and GS96 have suggested that it may occur at periods as short as 0^75 or O'^SO. Also, some 
clusters have bright variables with periods of about 0l5 and light curve morphology that is 
significantly different from RR Lyrae variables with the same period. Examples of such stars 
are V47, V68 and V123 in u Cen, V76 in M5 and VIO in NGC 5897. Nemec et al. (1994) 
have suggested that V68 in uj Cen might be an anomalous Cepheid or a first-overtone AHBl 
variable. There is also some uncertainty about classification at the longer periods. V2 in 
NGC 6712 has been classified as semi-regular by Rosino (1966, 1978) but Barnes & Dupuy 
(1975) considered it to be an RV Tauri variable. As high precision light curves become 
available for more stars, it may be possible to use light curve morphology, in addition to 
color and luminosity, to distinguish between the different types of pulsating variables. 

6. THE SR AND RED VARIABLES 

The SR and red variables are the pulsating giant stars that Rosino (1978) designated 
as yellow semi-regular, red semi-regular, irregular or Mir a variables in his review of variable 
stars in globular clusters. Whitelock (1986) made a plot of bolometric magnitude versus 
the log of the fundamental period for cluster variables with periods greater than 1 day and 
pointed out that the relationship for the yellow and red variables might form an extension 
to the population II Cepheid P-L relation. 

Since not all authors have designated a specific classification for variables with longer 
periods, we have included in column 7 of Table 1 the stars classified as semi-regular or 
irregular variables and all the stars with published periods greater than 35 days (excluding 
RV Tauri variables). These types of stars have not been well surveyed in the Galactic 
globular clusters because most search programs are designed to optimize detection and 
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period determination for variables with periods less than one day. Thus our sample is 
probably incomplete. Periods have been published for some, but not all of these stars 
because of the irregular nature of their variation. Many of the 'unclassified' variables 
for which periods have not been determined (column 3 of Table 1) may belong in this 
category. The period-frequency plot in the bottom panel of Figure 1 shows that only five 
globular cluster variables have periods greater than 300 days, but according to the General 
Catalogue of Variable Star^ (Kholopov 1985a,b, 1987), periods greater than 300 days are 
not uncommon among field stars. Also, in a study of Mira variables in the Galactic bulge, 
Whitelock (1990) reported periods ranging up to 720 days with most between 360 and 560 
days. This is probably a metallicity effect. Mira variables with such long periods do not 
occur in metal poor systems like globular clusters. However, in a study of the MACHO data 
for red giants in the LMC, Wood et al. (1999) have shown that there are red variables with 
periods up to about 1000 days. Because the MACHO observations of the LMC were made 
throughout the year, for a few years, variations on such a long time scale could be readily 
detected. Thus, as we have noted above, our sample of SR/red variables is undoubtedly 
incomplete. 
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FIGURE CAPTIONS 

Fig. 1. — Period- frequency distribution for pulsating variables in Galactic globular clusters. 
The upper panel shows the distribution for stars with periods less than 1 day. The stars with 
periods less than 0^20 are SX Phe variables. Most of the others are RR Lyrac variables, but 
some of the longer period variables may be anomalous Cepheids or population II Cepheids. 
In the middle panel, the Cepheids and RV Tauri variables are plotted. The periods plotted 
for the RV Tauri stars are their pulsation periods. Also included are two anomalous Cepheids 
(one in uj Cen and one in NGC 5466) with periods less than a day. In the lower panel, the 
distribution for pulsating variables with periods between 50 and 440 days is shown. 

Fig. 2. — Period-frequency distribution for eclipsing binary stars in the Galactic globular 
clusters. 

Fig. 3. — Period-frequency distribution for RR Lyrae variables in clusters of Gosterhoff type 
I (upper panel) and Gosterhoff type II (lower panel) . 

Fig. 4. — Plots of [Fe/H] (from Table 1) versus mean period for the globular cluster RRO 
variables. In the upper panel, only clusters with at least 15 RRO variables are included and 
in the lower panels clusters with at least 10 and at least 5, respectively are included. Filled 
circles represent clusters that have red horizontal branches, i.e. HBR=[(B-R)/(B-|-V-|-R)]< 
and open circles represent clusters with HBR> 0. Crosses represent clusters for which HBR 
has not been determined. (For ui Cen, a mean [Fe/H] value, —1.62, has been adopted.) 

Fig. 5. — Plots of HBR (from Table 1) versus mean period for the globular cluster RRO 
variables. In the upper panel, only clusters with at least 15 RRO variables are included and 
in the lower panels clusters with at least 10 and at least 5, respectively are included. Filled 
circles represent clusters that have [Fe/H]> —1.6 and open circles represent clusters with 
[Fe/H]< -1.6. 
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Fig. 6. — Plots of [Fe/H] (from Table 1) versus mean period for the globular cluster RRl 
variables. We include the RROl stars here as well because their dominant mode of pulsation 
is generally the first overtone. In the upper panel, only clusters with a total of at least 15 
RRl and RROl variables are included and in the lower panels clusters with at least 10 and 
at least 5, respectively are included. The symbols are the same as in Figure 4. 

Fig. 7. — Plots of HBR (from Table 1) versus mean period for the globular cluster RRl 
and RROl variables. In the upper panel, only clusters with a total of at least 15 RRl and 
RROl variables are included and in the lower panels clusters with at least 10 and at least 5, 
respectively are included. The symbols are the same as in Figure 5. 

Fig. 8. — Plots of [Fe/H] (from Table 1) versus mean fundamentalized period for the globular 
cluster RR Lyrae variables. In the upper panel, only clusters with at least 15 RR Lyrae 
variables are included and in the lower panels clusters with at least 10 and at least 5, 
respectively are included. The symbols are the same as in Figure 4. 

Fig. 9. — Plots of HBR (from Table 1) versus mean fundamentalized period for the globular 
cluster RR Lyrae variables. In the upper panel, only clusters withat least 15 RR Lyrae 
variables are included and in the lower panels clusters with at least 10 and at least 5, 
respectively are included. The symbols are the same as in Figure 5. 

Fig. 10. — Plots of [Fe/H] (upper panel) and HBR (lower panel) versus pulsation period for 
the individual Cepheids and RV Tauri variables in GGCs. In the upper panel, the symbols 
are the same as in Figure 4. In the lower panel, the symbols are the same as in Figure 5. 
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TABLE CAPTIONS 

Table 1. Numbers and Types of Variable Stars in Galactic Globular 
Clusters 

Table 2. Mean Periods of RR Lyrae Variables in Galactic Globular 
Clusters 
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Table 1. Numbers and Types of Variable Stars in Galactic Globular Clusters 



Cluster 


HBR/[Fe/H] 


No./ 


SX 


RR 


Cep or 


SR or 


Eel 


Ref 






(periods) 






RVTau 


Red 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 104/47TUC 


-0.99/-0.76 


50/(37) 


3 


1 





14 


19 


1,2,3,4 


NGC 288 


0.98/-1.24 


10/(10) 


6 


2 





1 


1 


5,6 


NGC 362 


-0.87/-1.16 


13/(8) 





7 





1 





7 


NGC 1261 


-0.71/-1.35 


20/(19) 





18 





1 





8,9 


NGC 1851 


-0.36/-1.22 


32/(29) 





29 











10,11,12,13 


NGC 1904/M79 


0.89/-1.57 


8/(3) 





3 





1 







NGC 2298 


0.93/-1.85 


4/(4) 





4 











14,15 


NGC 2419 


0.86/-2.12 


41/(33) 





31 


1 


4 





16 


NGC 2808 


-0.49/-1.15 


5/(3) 





2 


1 


2 





17 


Pal 3 


-0.50/-1.66 


12/(1) 








1 








18,19,20 


NGC 3201 


0.08/-1.58 


92/(77) 





77 











21,22 


Pal 4 


-1.00/-1.48 


2/(2) 











2 







NGC 4147 


0.55/-1.83 


16/(9) 





9 











23 


NGC 4372 


1.00/-2.09 


16/(13) 


8 





2 





3 


24 


Rup 106 


-0.82/-1.67 


16/(16) 


3 


13 











25 


NGC 4590/M68 


0.17/-2.06 


44/(44) 


2 


42 











26,27 


NGC 4833 


0.93/-1.79 


23/(15) 





14 





2 





28 


NGC 5024/M53 


0.81/-1.99 


67/(60) 





58 





2 





29,30 


NGC 5053 


0.52/-2.29 


15/(14) 


5 


9 











31 


NGC 5139/a; Gen 


— /-1.62 


275/(249) 


34 


161 


10 


15t 


29 


32,33,34 


















35,36,37 


NGC 5272/M3 


0.08/-1.57 


254/(187) 


1 


182 


1 


3 





38,39,40 



41,42,43 
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Table 1 — Continued 



Cluster 


HBR/[Fe/H] 


No./ 


SX 


RR 


Cep or 


SR or 


Eel 


Ref 






(periods) 






RVTau 


Red 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 5286 


0.80/-1.67 


24/(13) 





13 











44,45,46 


NGC 5466 


0.58/-2.22 


31/(30) 


6 


20 


1 





3 


47,48 


















49,50 


NGC 5634 


— /-1.82 


7/(6) 





6 











51 


IC 4499 


0.11/-1.60 


113/(98) 


1 


97 











52,53,54 


NGC 5824 


0.79/-1.85 


27/(7) 





7 













Pal 5 


-0.40/-1.43 


5/(5) 





5 













NGC 5897 


0.86/-1.80 


9/(9) 


1 


7 





1 





55,56,57 


NGC 5904/M5 


0.31/-1.29 


158/(135) 


5 


126 


2 


1 


1 


58,59,60,61 


















62,63,64 


















65,66,67 


NGC 5927 


-1.00/-0.37 


9/(5) 











5 





68 


NGC 5946 


— /-1.38 


6/(2) 





2 











69 


NGC 5986 


0.97/-1.58 


10/(8) 





8 





1 





70,71 


NGC 6093/M80 


0.93/-1.75 


8/(7) 





6 


1 








72 


NGC 6121/M4 


-0.06/-1.20 


67/(57) 


4 


40 





1 


8 


73,74,75,76 


















77,78,79 


NGC 6101 


0.84/-1.82 


15/(0) 

















80 


NGC 6144 


1.00/-1.73 


1/(0) 

















69 


NGC 6139 


0.91/-1.68 


10/(4) 





4 











81 


NGC 6171/M107 


-0.73/-1.04 


23/(22) 





22 











82,83 


NGC 6205/M13 


0.97/-1.54 


29/(18) 





5 


5 


8 





84,85,86 



87,87,89 
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Table 1 — Continued 



Cluster 


HBR/[Fe/H] 

/ L / J 


No./ 


SX 


RR 


Cep or 


SR or 


Eel 


Ref 






(periods) 






RVTau 


Red 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 6218/M12 


0.97/-1.48 


1/(1) 








1 








90,91 


NGC 6229 


0.24/-1.43 


47/(42) 





38 


2 





2 


92,93 


NGC 6235 


0.89/-1.40 


4/(4) 





3 








1 


94 


NGC 6254/MlO 


0.98/-1.52 


4/(2) 








2 


1 





95 


NGC 6266/M62 


0.32/-1.29 


89/(74) 





74 











96 


NGC 6273/M19 


— /-1. 68 


5/(5) 





1 


4 








97 


NGC 6284 


— /-1.32 


10/(9) 





6 


2 


2 





98 


NGC 6287 


0.98/-2.05 


3/(0) 

















99 


NGC 6293 


0.90/-1.92 


5/(4) 





4 





1 





100 


NGC 6304 


-1.00/-0.59 


21/(0) 

















101 


NGC 6341/M92 


0.91/-2.29 


21/(20) 


2 


17 


1 








102,103 


NGC 6333/M9 


0.87/-1.72 


20/(20) 





17 


1 


1 


1 


104 


NGC 6356 


-1.00/-0.50 


8/(3) 











3 





68,105 


NGC 6352 


-1.00/-0.70 


3/(1) 











1 





68,106,107, 


NGC 6366 


-0.97/-0.82 


1/(1) 





1 











108 


NGC 6362 


-0.58/-0.95 


47/(37) 


1 


35 








1 


109,110,111 


Haute PI 


— /-1. 50 


15/(0) 



















NGC 6380 


— /-0.50 


1/(0) 



















Terzan 1 


— /-0.35 


4/(0) 

















112 


NGC 6388 


— /-0.60 


29/(10) 





10 











68,113, 


















114,115 


Pismis 26 


— /-0.50 


2/(0) 



















NGC 6402/M14 


0.65/-1.39 


90/(61) 





54 


6 


6 





116 
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Table 1 — Continued 



Cluster 


HBR/[Fe/H] 

/ L / J 


No./ 


SX 


RR 


Cep or 


SR or 


Eel 


Ref 






(periods) 






RVTau 


Red 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 6401 


— /-I. 12 


25/(0) 

















112 


NGC 6397 


0.98/-1.95 


7/(6) 


2 











4 


117,118 


NGC 6426 


0.53/-2.26 


15/(14) 





14 











119,120 


Terzan 5 


-1.00/-0.28 


4/(4) 





1 





2 


1 


121,122 


NGC 6441 


— /-0.53 


37/(10) 





10 











115,123 


NGC 6522 


0.71/-1.44 


1/(1) 





1 











124 


NGC 6539 


-1.00/-0.66 


1/(0) 



















NGC 6544 


1.00/-1.56 


2/(1) 





1 











125 


NGC 6541 


1.00/-1.83 


1/(0) 

















126 


NGC 6553 


-0.34/- 1.00 


13/(3) 





2 





1 







NGC 6558 


— /-1. 44 


15/(9) 





9 











127,128 


IC 1276/Pal 7 


— /-0.73 


5/(1) 





1 





4 







NGC 6569 


~ /-0.86 


19/(0) 

















129 


NGC 6584 


-0.15/-1.49 


46/(42) 





42 











130,131 


NGC 6624 


-1.00/-0.42 


4/(0) 

















132,133 


NGC 6626/M28 


0.90/-1.45 


20/(15) 





10 


3 


4 


2 


134,135,136 


NGC 6638 


-0.30/-0.99 


61/(19) 





13 





6 





137 


NGC 6642 


— /-1.35 


18/(16) 





16 











138 


NGC 6637/M69 


-1.00/-0.71 


7/(2) 











2 







NGC 6652 


-1.00/-0.96 


5/(0) 

















139,140 


NGC 6656/M22 


0.91/-1.64 


36/(24) 





18 


1 


5 





141,142,143 


NGC 6681/M70 


0.96/-1.51 


3/(2) 





2 











144,145 


NGC 6712 


-0.64/-1.01 


19/(14) 





9 





5 





146 
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Table 1 — Continued 



Cluster 


HBR/[Fe/H] 

/ L / J 


No./ 


SX 


RR 


Cep or 


SR or 


Eel 


Ref 






(periods) 






RVTau 


Red 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 6715/M54 


0.87/-1.59 


94/(64) 





61 


2 





1 


147 


NGC 6717 


0.98/-1.29 


1/(1) 





1 











148 


NGC 6723 


-0.08/-1.12 


32/(28) 





28 





1 





68,149,150 


NGC 6752 


1.00/-1.56 


10/(7) 


3 





1 





3 


151,152 


















153,154 


NGC 6760 


-1.00/-0.52 


4/(0) 



















NGC 6779/M56 


0.98/-1.94 


7/(4) 





2 


2 








155 


NGC 6809/M55 


0.87/-1.81 


37/(37) 


24 


13 











156,157 


NGC 6838/M71 


-1.00/-0.73 


26/(18) 


3 








2 


13 


88,158 


NGC 6864/M75 


-0.42/-1.32 


17/(6) 





6 





3 





159 


NGC 6934 


0.25/-1.54 


86/(82) 


1 


78 





1 


2 


160,161 


NGC 6981/M72 


0.14/-1.40 


35/(25) 





25 











162 


NGC 7006 


-0.28/-1.63 


72/(64) 





62 





1 





163,164 


NGC 7078/M15 


0.67/-2.25 


154/(104) 


2 


88 


2 





1 


165,166,167 


















168,169,170 


















171,172 


NGC 7089/M2 


0.96/-1.62 


33/(33) 





29 


4 








173 


NGC 7099/M30 


0.89/-2.12 


13/(3) 





3 











112 


Pal 12 


-1.00/-0.94 


3/(0) 



















Pal 13 


-0.20/-1.65 


4/(4) 





4 













NGC 7492 


0.81/-1.51 


4/(4) 





3 


1 










TOTAL 




2993(2220) 


117 


1842 


60 


117 


96 





-39- 



fo; Cen contains 5 spotted variables with periods that range from 3.3 days to 34 days. 
It also contains 2 possible ellipsoidal variables. Since variables of these types have only been 
identified in uj Cen, they are included in this column, but they are not plotted in Figure 1. 
There are only 3 SR/red variables in u Cen. 
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Table 2. Mean Periods of RR Lyrae Variables in Galactic Globular Clusters 





#of 


<P> 


#of 


<P> 


#of 


<P> 


#of 


<P> 


<Pf> 


Cluster 


RRO 


(RRO) 


RRl 


(RRl) 


RR2 


(RR2) 


RROl 


(RROl) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 104 


1 


0.737 

















0.737 


NGC 288 


1 


0.678 


1 


0.430 












0.628 


NGC 362 


7 


0.542 

















0.542 


NGC 1261 


13 


0.555 


5 


0.328 












0.523 


NGC 1851 


21 


0.571 


7 


0.317 


1 


0.266 







0.531 


NGC 1904 


2 


0.685 


1 


0.335 












0.607 


NGC 2298 


1 


0.640 


3 


0.394 












0.557 


NGC 2419 


24 


0.655 


6 


0.380 







1 


0.407 


0.624 


NGC 2808 


1 


0.539 







1 


0.306 







0.526 


NGC 3201 


72 


0.554 


3 


0.338 


2 


0.274 







0.547 


NGC 4147 


4 


0.531 


5 


0.339 












0.489 


Rup 106 


13 


0.617 

















0.617 


M68 


13 


0.613 


17 


0.390 







12 


0.397 


0.554 


NGC 4833 


7 


0.708 


7 


0.380 












0.609 


M53 


29 


0.649 


27 


0.351 


2 


0.309 







0.562 


NGC 5053 


5 


0.672 


4 


0.354 












0.585 


u! Cen 


76 


0.651 


59 


0.387 


26 


0.306 







0.581 


M3 


145 


0.555 


27 


0.345 


5 


0.286 


5 


0.358 


0.537 


NGC 5286 


8 


0.614 


5 


0.344 












0.556 


NGC 5466 


13 


0.646 


5 


0.382 


2 


0.262 







0.592 


NGC 5634 


3 


0.621 


3 


0.379 












0.565 


IC 4499 


63 


0.580 


14 


0.354 


3 


0.284 


17 


0.358 


0.544 


NGC 5824 


7 


0.624 

















0.624 
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Table 2 — Continued 





# of 


<P> 


# of 


<P> 


# of 


<P> 


# of 

TT 


<p> 


<Pf> 


Cluster 


RRO 


(RRO) 


RRl 


(RRl) 


RR2 


(RR2) 


RROl 


(RROl) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Pal 5 







3 


0.316 


2 


0.273 







0.438 


NGC 5897 


3 


0.828 


2 


0.437 


2 


0.345 







0.688 


M5 


91 


0.551 


34 


0.320 


1 


0.265 







0.517 


NGC 5946 


2 


0.651 

















0.651 


NGC 5986 


7 


0.652 


1 


0.328 












0.626 


M80 


4 


0.651 


2 


0.366 












0.598 


M4 


31 


0.533 


7 


0.302 


2 


0.250 







0.505 


NGC 6139 


3 


0.662 


1 


0.417 












0.636 


M107 


15 


0.538 


7 


0.287 












0.490 


M13 


1 


0.750 


3 


0.358 


1 


0.313 







0.543 


NGC 6229 


30 


0.553 


7 


0.330 


1 


0.264 







0.529 


NGC 6235 


2 


0.600 


1 


0.352 












0.557 


M62 


62 


0.544 


11 


0.304 


1 


0.248 







0.522 


M19 


1 


0.507 

















0.507 


NGC 6284 


6 


0.588 

















0.588 


NGC 6293 


2 


0.600 


2 


0.353 












0.537 


Mq2 


1 1 

J. J. 


630 


(J 


364 


1 

J. 


31 3 


n 

yj 




582 


M9 


8 


0.638 


9 


0.342 












0.543 


NGC 6366 


1 


0.513 

















0.513 


NGC 6362 


18 


0.547 


15 


0.299 


2 


0.251 







0.477 


NGC 6388 


5 


0.712 


5 


0.332 












0.579 


M14 


39 


0.564 


10 


0.335 


5 


0.320 







0.540 


NGC 6426 


9 


0.704 


4 


0.332 







1 


0.408 


0.619 
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Table 2 — Continued 





# of 


<P> 


# of 


<P> 


# of 


<p> 


# of 

TT 


<p> 


<pf > 


Cluster 


RRO 


(RRO) 


RRl 


(RRl) 


RR2 


(RR2) 


RROl 


(RROl) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Terzan 5 


1 


0.60 

















0.60 


NGC 6441 


10 


0.768 

















0.768 


NGC 6522 


1 


0.564 

















0.564 


NGC 6544 


1 


0.57 

















0.570 


NGC 6553 


2 


0.526 

















0.526 


NGC 6558 


6 


0.556 


3 


0.345 












0.525 


Pal 7 


1 


0.548 

















0.548 


NGC 6584 


34 


0.560 


8 


0.304 












0.531 


M28 


8 


0.577 


2 


0.312 












0.545 


NGC 6638 


1 


0.666 


11 


0.307 


1 


0.308 







0.439 


NGC 6642 


10 


0.544 


6 


0.322 












0.502 


M22 


10 


0.632 


8 


0.361 












0.567 


M70 


1 


0.564 


1 


0.402 












0.552 


NGC 6712 


7 


0.557 


2 


0.338 












0.534 


M54 


55 


0.579 


6 


0.342 












0.568 


NGC 6717 


1 


0.575 

















0.575 


NGC 6723 


23 


541 


4 


292 


1 

J. 


288 


n 




51 7 


M56 


1 


0.906 


1 


0.423 












0.736 


M55 


4 


0.662 


5 


0.391 


4 


0.321 







0.571 


M75 


3 


0.531 


3 


0.271 












0.447 


NGC 6934 


68 


0.574 


9 


0.308 


1 


0.247 







0.553 


M72 


24 


0.547 


1 


0.353 












0.544 


NGC 7006 


53 


0.569 


9 


0.329 












0.550 
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Table 2 — Continued 





#of 


<P> 


#of 


<P> 


#of 


<P> 


#of 


<P> 


<Pf> 


Cluster 


RRO 


(RRO) 


RRl 


(RRl) 


RR2 


(RR2) 


RROl 


(RROl) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


M15 


39 


0.637 


30 


0.337 


2 


0.290 


17 


0.401 


0.551 


M2 


17 


0.725 


9 


0.345 


3 


0.298 







0.621 


M30 


3 


0.698 

















0.698 


Pal 13 


4 


0.572 

















0.572 


NGC 7492 


1 


0.805 


1 


0.292 


1 


0.280 







0.556 


TOTAL 


1269 


0.585 


447 


0.345 


73 


0.296 


53 


0.383 


0.550 
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